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Abstract 

We present a detailed discussion of a novel dynamical renormaliza- 
tion group scheme: the Dynamically Driven Renormalization Group 
(DDRG). This is a general renormalization method developed for 
dynamical systems with non-equilibrium critical steady-state. The 
method is based on a real space renormalization scheme driven by 
a dynamical steady-state condition which acts as a feedback on the 
transformation equations. This approach has been applied to open 
non-linear systems such as self-organized critical phenomena, and it 
allows the analytical evaluation of scaling dimensions and critical ex- 
ponents. Equilibrium models at the critical point can also be consid- 
ered. The explicit application to some models and the corresponding 
results are discussed. 



Key words: Renormalization group; non-equilibrium steady states; driven 
dynamical systems; self-organized criticality. 
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1 Introduction 

The study of second order phase transitions introduced in the field of sta- 
tistical physics the concept of criticality 0, H, H, El- A critical system is 
characterized by the absence of a characteristic lengthscale: the system fluc- 
tuates strongly and correlation functions show non-analytic behavior. In the 
critical region therefore the usual tools of a physicist, principally perturbation 
theory, fail completely 0. The renormalization group (RG) theory provides 
a comprehensive understanding of second order phase transitions and critical 
phenomena || || [7| ||. This theory has been a major breakthrough in statis- 
tical physics and lead to the study of many others scale invariant and critical 
phenomena. Among these phase transitions associated with non-equilibrium 
states H [UJ, fractal growth phenomena [PI , self-organized critical (SOC) 



systems |T2|, O] and a vast class of complex systems in which information 



spread over a wide range of length and time scales flU 



Since both areas of critical and complex systems deal with self-similar 
structures, it was natural for physicists who were familiar with RG techniques 
to consider these new problems as a possible playground for these methods. 
The situation turned out to be more complex, since the properties of non- 
equilibrium critical phenomena are quite different from those of ordinary 
critical systems. For instance, there is no ergodic principle and in general 
it is not possible to assign a Boltzmann weight to a configuration. For the 
above reasons many authors concluded that these problems pose questions 
of new type for which it would be desirable to have a common theoretical 
framework In the past decade an intense activity has been devoted to 
a better understanding of these systems and theoretical methods are being 



developed Ml 16 



Recently we introduced a renormalization scheme [O] for sandpile models 



[ p~3|1 , that has later been applied |L8] to forest-fire models |19], |2T)fl . This ap- 
proach deals with the critical properties of the system by introducing in the 
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renormalization equations a dynamical steady state condition which provides 
the non-equilibrium stationary statistical weights to be used in the calcula- 
tion. In this way it is possible to characterize the fixed point dynamics and 
to compute analytically the critical exponents. 

Here we present the general formulation of this novel type of dynami- 
cal renormalization group: the Dynamically Driven Renormalization Group 
(DDRG), a general theoretical method for dynamical non-equilibrium sys- 
tems with critical stationary state. The essential idea is to combine a real 
space renormalization group (RSRG) scheme with the dynamical steady-state 
condition which characterizes the stationary regime. The RG equations are 
driven by the steady state condition feedback from which we obtain the 
configurations approximate statistical weight to be used in the dynamical 
renormalization of the master equation. While the approximate station- 
ary distribution neglects correlations, these are considered in the dynamical 
renormalization. Finer scale correlations are then included in the approx- 
imate stationary distribution at the new scale which is calculated through 
the steady state condition with renormalized dynamical parameters. This 
strategy gives an RG scheme which can be improved considering increas- 
ingly better approximations in both the dynamical renormalization and the 
stationary distribution. 

For the sake of clarity we discuss extensively the explicit application of 
the method to some driven dissipative systems referring to the present gen- 
eral framework. In particular we report in full details the results obtained 
for the critical height sandpile model and the forest fire model. For this class 
of models the DDRG can be considered a general renormalization scheme 
which provides a new class of analytical tools for the study of the stationary 
critical state. The DDRG can also be applied to ordinary dynamical critical 
phenomena for which the stationary state is characterized by Gibbs distri- 
butions. In this case the DDRG can be directly compared with other RG 
approaches. 
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The paper is organized as follows: sec. 2 introduces the concept of non- 
equilibrium steady-state and the approximate description in term of mean- 
field stationary conditions. Sec. 3 presents in full detail the Dynamically 
Driven Renormalization Group. The basic recursion relations are obtained 
and the conceptual scheme is discussed as well as the approximations in- 
volved. Sec. 4 shows the explicit application of the DDRG to sandpile and 
forest fire models. Sec. 5 shows the application scheme to ordinary critical 
phenomena. Sec. 6 is devoted to conclusions and perspectives. In the ap- 
pendix A we discuss in detail the coarse graining of time. 

2 Equilibrium and non-equilibrium steady states 

The distribution over configuration space in equilibrium ensembles as well 
as many extended dissipative and nonlinear dynamical systems evolves in 
stationary states described by time independent probability distributions. 
The stationary state in some cases shows long range correlations and self- 
similar properties^]. To describe these phenomena several models have been 
introduced: here we will consider models defined on a discrete <i-dimensional 
lattice. To each site of the lattice is associated a variable Oi that can assume q 
different values (<7j = 1, 2, 3, - • • , g). The subscript i labels the lattice site. A 
dynamics characterized by a set of parameters^ \i acts on these variables and 
defines the temporal evolution of the model. The system can be described 
by P(<x, £), the probability that at the instant t the system is in the state 
o = {o"j}. This is the usual way to characterize physical ensembles for which 
we want to know the statistical distribution in phase space. 

For stationary processes, the system is invariant under uniform time 
translations. Thus, the variables <Ji(t + r) are statistically indistinguish- 

^dn equilibrium ensembles this happens for special values of the control parameter, i.e., 
at critical points 

2 We should use a vector notation /i.= {/j,} to denote the set of all variables. For the 
sake of simplicity we use the simpler notation \x. 
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able from the untranslated variables (Ji{t). This implies that all single time 
averages are constants and there is a single time probability density W{o) 
independent of time. An ensemble of systems at thermodynamic equilibrium 
is, for example, stationary as well as an ensemble in which the component 
systems are maintained in a non-equilibrium steady-state, at least after the 
ensemble has been sufficiently aged. For equilibrium ensembles, the steady 
state statistical distribution is given by the Gibbs distribution. In this case 
the functional dependence of the equilibrium distribution function W(a) on 
the parameters^ fi should be consistent with statistical mechanics and prop- 
erly describe the static equilibrium properties of the system. In particular, 
the partition function and the equilibrium expectation values should show 
critical point singularities appropriate to the spatial dimensionality and the 
symmetry of the order parameter. For non-equilibrium ensembles in prin- 
ciple one does not know how to assign a statistical weight to a given non- 
equilibrium configuration. An additional problem is that, since the system 
is not in equilibrium, the distribution changes in time. 

One possibility is to study systems which have "settled down" into non- 
equilibrium steady states, so that the distributions, while non Gibbsian, have 
become stationary. In this case we can describe, at least approximately, the 
single-time statistical distribution |2l| , since the densities p K of sites in a 
particular state <Ji — k (k — 1, • • - , q) do not change, on the average, as a 
function of time. Associated to the set of time independent average densities 
{(Pk)} there is a unique stationary probability distribution, that character- 
izes single-time averages in the steady-state ensemble. We can therefore 
describe the average statistical state by means of stationarity conditions for 
the system. These can be obtained from mean field equations of the form 



where the operator describes the evolution of the system as a function of 
3 For instance the temperature or the set of applied fields. 



d_ 

Of 



(1) 
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the dynamical parameters defined above. For example, in the case of dilute 
gases the non-linear function that gives the rate of change of the density 
becomes a non-linear functional, i.e. the collision operator. In general, the 
operator iS M represents the sum of dissipative contributions, mechanical (non- 
dissipative) contribution, and the effect of other external fluxes. Note that 
the operator on the right-hand side of Eq. (jH) is time independent, so that 
the differential equation is first order in time. 

Time independent solutions of Eq. (H) will be referred to as "steady- 
states", although we should keep in mind that those are only the average 
states of the ensemble. Mathematically the steady-states (p K ) are determined 
by the equation 

^{{Pn)} = S,({( Pk )}) = (2) 

The above equation can have more than one solution, even when «S /Lt is a linear 
function. In the following, however, we will consider only the presence of a 
single meaningful stationary state. In ordinary statistical systems, the Eq. |2] 
represents the thermodynamic equilibrium condition. For driven dynamical 
systems, it describes the driving of the system to the non-equilibrium steady 
state, by means of a balance condition. We can express the stability or 
self-organization properties in mathematical form as 

fm(p K (t)) = (p K ) (3) 

for any initial condition of the Eq.([I]): the steady-state is an attractor for 
the dynamics. 

It is possible to go further on these lines to find a more accurate de- 
scription of the steady-state. For example, one could write down mean-field 
equations as Eq.(^j) for the average density correlations. We are, however, 
interested in systems which show critical properties in the stationary state. 
This implies non-analyticity and long-range correlations, so that it is impos- 
sible to go much further along with a mean-field description. To describe 
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critical systems and scaling behavior, we have to turn our attention to renor- 
malization group methods. 



3 The Dynamically Driven Renormalization 
Group 

In this section we present the general formalism of a new type of dynami- 
cal renormalization group especially suited for systems with non-equilibrium 
critical steady-state. In fact, by using the dynamical steady-state condition 
we are able to develop a renormalization strategy which allows us to compute 
critical exponents in a wide range of non-equilibrium systems. The method 
can also be applied to equilibrium dynamical critical phenomena for which 
the steady-state distribution is given by the equilibrium Gibbs distribution. 

In the following we will refer to RSRG schemes. The real-space formula- 
tion of the renormalization group, by virtue of its simplicity and versatility 
is a vital tool in the theoretical understanding of critical phenomena. In 
contrast to momentum-space renormalization group, the RSRG, in general 
does not have any systematic way to treat the approximations: there is no 
small parameter that controls an expansion. On the other hand, to improve 
the accuracy of the method one can use "higher order" techniques like the 
introduction of proliferation (additional couplings), larger cells, toroidal or 
rectangular transformation or other extrapolation scheme 0, [g, In this 
way, the critical exponents for several models have been obtained with good 
precision. 



3.1 Coarse graining and renormalization 

The essential ideas of the dynamic real space renormalization group approach 
derive from Kadanoff's block analysis |23| and from coarse graining of time 
proposed by Suzuki |24]. We begin summarizing the main derivation of the 
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method. We start with the following general Master Equation (ME): 

^P(a,t) = jC(n)P(a,t) (4) 

where P(a, t) denotes the probability distribution function for the configura- 
tions of the system at time t, and £(fi) is the temporal evolution operator, 
characterized by a set of dynamical parameters /i. With the operator TZ 
we indicate the coarse graining operator, that eliminates degrees of freedom 
inside a cell and rescales time and space. The application of 1Z yields 

KP(a,t) = P'(S,t') (5) 

where P'(S,t') denotes the probability distribution for the coarse grained 
system. More explicitly, we can write: 

KP(a, t) = n(e tc P(a, 0)) = e t,c 'P'(S, 0) (6) 



The scale transformation tC — > t'C yields the dynamical RG approach 



while the scale transformation of P(a, 0) corresponds to the usual static RG 
approach [H, [7]]. Denoting with the vector fi the parameters of the system, 
the RG yields recursion relations: 

// = /(//) and t' = g{t, //) ~ t g(n) (7) 

from which it is possible to calculate the fixed points and the critical expo- 
nents of the model. These equations are obtained from the renormalization 
procedure which impose that P'(S, t') has the same functional form as P(cr, t). 

The most delicate problem in RSRG approaches is to take a partial elim- 
ination of degrees of freedom in Eq.s (H) and (H) f|- To deal with this 



problem many approximate methods have been proposed p5|, |26[ and have 
been mostly applied to the kinetic Ising model ^7| . 



Achaim and Kosterlitz |28[ perform a Migdal-Kadanoff |29[ transforma- 
tion on the probability distribution function, assuming a functional form of 
the type: 

p^ t ) = e -H(«)+h(t) (g) 
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and treating the time dependent part as a perturbation. In this way they 
study the relaxation of the probability density towards equilibrium. The 
dynamical critical exponent is extracted from the scaling of the slower relax- 
ation mode. Suzuki et al.|3(| obtain the coarse graining of time studying the 



dynamical equation for the first moment of the probability distribution (i.e. 
the magnetization). The equation is decimated in the Migdal-Kadanoff ap- 
proximation and the equation of motion for the decimated modes are solved 
perturbatively. The time scaling is chosen so that it keeps the equations 
in the same form. A different formalism, suitable to study the properties 
of the model even far from the critical region, was developed by Mazenko 
et al.|31| The coarse graining operator and the time rescaling are chosen 
self-consistently in order to insure the Markoffian behavior of the renormal- 
ized spin flip operator. One then writes recursion relations for the two-point 
correlation functions from which the critical properties of the model are ex- 
tracted. 

To develop a RSRG method suitable for irreversible non-equilibrium sys- 
tems we will consider a more explicit treatment of the master equation. In 
particular we focus our attention on the dynamics of discrete models on a 
lattice characterized by a set of lattice variables a = {o"i}, each of them being 
a g-state variable (see previous section), and by the dynamical parameters 
fi. In this case we can rewrite eq.(|j) as 

P(a,t) = J2(v\T(ri\a )P(a°,0) (9) 

where (a | T(/i) | a ) is the transition probability from the configuration 
a = {er°} to the configuration a = {o~i} in a unit time step t. The symbol 
Y^{a} wm always mean a summation over all the configurations. The operator 
T is the discrete counterpart of the operator C 

We then coarse grain the system by rescaling lengths and time according 
to the transformation x — > bx and t — > b z t. The renormalization transfer- 
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mation can be constructed through a renormalization operator 1Z(S, er) that 
introduces the coarse grained variables set S = {Si} and rescales the lengths 
of the system This operator contains all the information connecting the 
coarse-grained state with the original one. Not every choice of the operator 
TZ will lead to a meaningful transformation, and in constructing it one should 
be guided by physical insight. Moreover, the transformation should satisfy 
some general properties in order to respect the dimension and the symmetry 
of the internal space of the dynamical variables, i.e. the renormalized vari- 
ables should be of the same kind as the original. In general, TZ is a projection 
operator with the properties: 

K(S, a) > for any S = {Si}, a = {a,} (10) 

and 

J2nS,a) = l. (11) 
{s} 

These properties preserve the normalization condition of the renormalized 
distribution. The explicit form of the operator TZ will be defined case by case 
in the various application of the method. Usually, it corresponds to a block 
transformation in which lattice sites are grouped together in a super-site that 
defines the renormalized variables S{ by means of a majority or spanning rule. 

We subdivide the time step in intervals of the unitary time scale and we 
apply repeatedly the dynamical operator T, obtaining 

P{a,N) = y E{<T\T N (n)\o o )P{o°,0) (12) 

{<T°} 

where T N denotes the application of the T operator N times. We can there- 
fore write the eq(0) for the coarse graining of the system as follows: 

p'(s, o = £ ns, <y) E <* I TbZ (A*) I °°)P(v°, 0) (13) 
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where we have included the application of the operator TZ and t' = b z t. The 
meaning of (cr \ T v (fi) | a ) has to be defined explicitly. In the simplest 
case h z = N where N is an integer number. Whenever is possible to define a 
continuous time evolution operator, T feZ can be found as an integration over 
time. In general, since we are often dealing with discrete time evolution, 
we have to consider the T bz as an effective dynamical operator. We will 
specify a projection operator T> for the dynamics, that samples only the paths 
which lead to an appropriate definition of the dynamical process at the coarse 
grained scale. Also in this CclSC, clS for the operator 1Z, spanning conditions 
or majority rules are obtained from physical considerations. The projection 
operator is chosen in such a way to preserve the form of the operator T at 
every scale. In this way it is possible to define recursion relations for the 
dynamical parameters. In appendix A, we present in detail the definition 
of the dynamical projector operator and the explicit form of the effective 
operator T bz . 

In order to define the RG transformation more clearly, eq. ( |T3| ) can be 
written as: 

P'(S,t') = E E n(S°,a°)J2n(S,a){a I T fcZ (/i) I a°)P(a°,0) (14) 

{a°}{S°} W 

where we used the properties X^s } 7l(S , a ) = 1. By multiplying and 
dividing each term by P'(S°, 0) = 7l(S°, a°)P(a°, 0), and changing the 

order of summations we have 

P'(S,t') = 

^{s°}{ £ {ff0} rc(sV°)P( < T°,o) J^,U) {lb) 

which finally identifies the renormalized dynamical operator T' and yields 
the coarse grained master equation in the form 

p'(S, t') = J2( S \ t '(a0 I S°)P'(S°, 0) (16) 

{SO} 
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In the following, we apply this scheme to systems with a stationary dis- 
tribution P(<J,t — > oo) = W{o) given either by an equilibrium or a non- 
equilibrium steady-state. We can therefore study the renormalization of the 
dynamics for small deviations from the stationary state 

P(a,t) = W(a) +eh(a,t) (17) 

where W(p) is the steady-state probability distribution and e is an expansion 
parameter for the non-stationary part of the distribution. In the lowest order 
approximation we can consider only the stationary part of the distribution 
in the renormalization of the evolution operator T, which, by comparing 
Eq. (|15|) and fllED, is given by: 

(S | T'(pt) | S°) = 

Z W o } J2 W} n(S ,a )K(S,a)(a \ T b > | <t°)W{o°) 

£{<x°} n{S°,a°)W{a°) • 1 ' 

This is the basic renormalization equation that defines the dynamical evo- 
lution operator for the coarse grained system. In principle one could also 
consider higher orders in e with an 74nsai20 for the form of h. For the sys- 
tems under considerations the zero order approximation will give already non 
trivial results, since we are interested in the dynamics of the system once the 
steady state has been reached. 

To understand intuitively the above transformation (Eq.[L8]) we can con- 
sider that the operator 7Z(S, a) is a projection operator that selects only the 
4 To first order in e one obtains: 

h'(S, = ]T ((S | T'(n) I S°) h'(S°, 0) + I T'( M ) | S°) E=Q W'(S°) 
{s°} £ 

where the dependence of T on e has been made explicit. This equation describes the 
relaxation to the stationary state from a non stationary configuration. Compare this with 
p8| | where the relaxation to the equilibrium state after a perturbation (i.e. magnetic field) 
was studied in a similar way. 
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configurations {&{} which, at the coarse grained level, are mapped into the 
configuration {Si}. The right-hand side term of the above equation can be 
read as follows. The operator TZ(S°, a ) selects only the configurations which 
are renormalized in the starting configuration {Sf}, each of them with a rel- 
ative weight given by W{o°) normalized with the factor Y,^(S°,a )W(a°). 
For each of these configurations we compute the statistical weight of the evo- 
lution to the configuration through the paths selected by the dynamical 
operator T v (see the appendix). Finally, we sum up only the contribution of 
the configurations {a{\ that renormalize in the configuration {Si}, selected 
from the operator 1Z(S, a). In other words, the new dynamical operator T' is 
defined as the sum of the statistical weight of all the paths that lead from any 
starting configuration that renormalizes in {Sf}, to any final configuration 
that renormalizes in {Si}. Each path has then an additional weight given by 
the relative probability of the starting configuration {of}. This last factor is 
the stationary statistical distribution, and preserve the normalization of the 
transition probability matrix T' . It is interesting to note that the stationary 
configurations distribution compares explicitly as a statistical weight in the 
renormalization equations. 

The Eq. (|TE|) is the basic renormalization equation from which the desired 
recursion relations are obtained. Imposing that the renormalized operator 
T' has the same form of the operator T, i.e. T'(fi) = T(fi'), we obtain 
the rescaled parameter set fi' = /(/i). This implies that the renormalized 
single time distribution P'(S, t') has the same functional form of the original 
distribution P(a, t). Since we are dealing with discrete evolution operators T, 
we define the time scaling factor b z as the average number of steps we apply 
the operator T in order to obtain that T'(/i) = T(//) for the coarse grained 
system. It therefore depend upon the parameters set \i. This condition 
defines the time recursion relation g(fi), from which it is possible to calculate 
the dynamical critical exponent z (see Eq. |7|). 
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3.2 Driving condition and recursion relations 

The scheme discussed so far is a general formulation valid for each system 
which exhibits a stationary state, and its application presupposes the knowl- 
edge of the explicit form of the steady-state distribution W. For example, in 
equilibrium phenomena W is given by the Gibbs distribution. In that case 
it is possible to apply several methods such as cumulant expansions and ex- 
act or approximate decimation to obtain the form of the recursion relations. 
For non-equilibrium dynamical systems in general we do not know the form 
of the steady-state distribution. We will therefore develop an approximate 
method to evaluate the stationary distribution to be used in the calculation 
of the renormalized master equation. 

The steady-state distribution can in general be split in two parts 

W(a) = WW (a) + (a) (19) 

where W^(a) and (er) are, respectively, the incoherent and coherent 

part of the distribution. The incoherent part of the distribution has the 
property 

Y,°i<TjW { i i X°) = &- (2°) 

M 

where a is the single site average. Hence, it does not include correlations 
among variables and expresses a mean field approximation for the system. 
The coherent part w( c \a) can be subdivided in parts describing different 
kind of correlations: nearest-neighbors, next-nearest-neighbors etc. The in- 
coherent part is a factorized distribution, that, for systems characterized by 
a q-state variables (see sect. 2), has the form 

W^(a) = U(p at ) (21) 

i 

where (p K ) is the average density of sites in the At-state. In this way, we 
have approximated the probability of each configuration {<7j} as the product 
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measure of the mean field probability to have a state <Tj in each corresponding 
site. The incoherent part contribution to the renormalization equation is 
therefore particularly easy to obtain: for non-equilibrium steady-state we 
can use the stationarity condition, S^({(p K )}) = 0, to evaluate the densities 
(p K ). By solving the stationary condition equation, the average densities of 
the K-states for the coarse grained system are obtained as a function of the 
dynamical parameters p at the corresponding iteration of the RG equations. 
By inserting this approximate distribution in Eq. (|18D , we get the following 
set of renormalization equations 



(S | T(p') | S°) 

4<j0} J2{a} 



(22) 



EK } ft(S*\<7°)IL<Po?> 

S^Opk)}) = (23) 

where the second equation denotes the dynamical steady state condition 
that allows the evaluation of the approximate stationary distribution at each 



coarse graining scale. We call the Eq. fl23|) the driving condition, since it 
drives the RG equations acting as a feedback on the scale transformation. 

Rewriting these equations in the form of recursion relations and adding 
the equation for the rescaling of time (see appendix A) we obtain (]: 

/*' = /(/*, {(ft.)}) ( 24 ) 
t' = tg(p,{(p K )}) (25) 
SMp K )})=0 (p K ) = u K (ji) (26) 

in which the driving condition appears explicitly. The above set of equa- 
tions (or, equivalently, the Eq.s (|22"|)-(|2"3|) synthesizes the Dynamically Driven 
Renormalization Group at the lowest order. The fixed points p* = f{p*) of 

5 Note that the /x recursion relation has to be read as a vectorial equation for the 
complete dynamical parameters set. 
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the first equation govern the critical behavior of the system. The second 
equation gives the dynamical critical exponent z: 



The third equation, the driving condition, can be seen as a feedback mech- 
anism between the dynamical and stationary properties of the system. In 
addition, from the fixed point solutions 



the average stationary properties for the critical state are obtained. 

In this form of the DDRG, we take into account only the uncorrelated 
part of the steady-state probability distribution. The results obtained are 
non trivial because correlations in the systems are considered in the dynam- 
ical renormalization of the operator T, that given a starting configuration 
traces all the possible paths leading to the renormalized final configuration. 
Moreover, geometrical correlations are treated by the operator 1Z that maps 
the system by means of spanning conditions or majority rules. The renor- 
malized uncorrelated part of the stationary distribution is evaluated from 
the stationary condition with renormalized parameters, thus providing an 
effective treatment of correlations. 

In principle we can also refine the method by including higher order con- 
tributions to the unknown steady-state distribution. We have considered the 
simplest approximation for the driving condition, that takes into account only 
the single point densities. It is possible to consider also two points occupation 
probabilities, or even higher order clusters. One would then write general- 
ized mean field equations [[3^] for the n— point probability distributions and 
couple their solution to the renormalization of the dynamical operator. In 
other words, at every iteration of the scale transformation one should solve 
the generalized mean-field equations with renormalized couplings //. This is 




(27) 



z = 



log b 




(28) 
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in the spirit of RSRG calculations where clusters of spins of higher order are 
introduced progressively in the calculations . 

The DDRG scheme can be applied to both equilibrium and non-equilibrium 
critical systems, but it is particularly suitable to the latter case. In fact, the 
renormalization procedure does not act directly on the stationary probabil- 
ity distribution, which is in general unknown in non-equilibrium phenom- 
ena. The renormalization equations depend on the stationary probability 
distribution only to weight local configurations, allowing the use of vari- 
ous approximations for its evaluation. The relevant difference between our 
scheme and other dynamical RG methods is the fact that we obtain a set 
of equations which are independent on the specific form of the stationary 
probability distribution of the system. This perspective is quite different 
from several previous real space dynamical RG approaches which were based 
on the explicit knowledge of the stationary distribution or the detailed bal- 



ance hypothesis [31]. The application range of these methods was therefore 
restricted to to the relaxation dynamics of equilibrium systems. 

Our goal is to describe far from equilibrium critical systems for which the 
Gibbsian equilibrium description is not valid. The method therefore finds 
potential application to the wide range of non-equilibrium critical phenom- 
ena, such as driven-diffusive systems, cellular automata and contact processes 
[ |10|] . A real space renormalization treatment of these models appears to be 
particularly suitable since they are usually defined on a lattice with discrete 
time steps. On the other hand, the method is limited by the possibility to 
describe the system in terms of a reasonable mean field theory. In addition 
the present real space formulation is not easily implemented in systems that 
show non-local interactions. 

In the next section we present the explicit implementation of the DDRG 
to two specific cases. The purpose of the following section is to provide an 
example of how the method works in practice for non-equilibrium systems. 
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4 DDRG applied to sandpile and forest-fire 
models 

Many extended dissipative dynamical systems form structures with long 
range spatial and temporal correlations. The concept of self-organized criti- 
cality (SOC) has been invoked by Bak, Tang and Wiesenfeld |13| to describe 
the tendency of slowly driven systems to evolve spontaneously toward a crit- 
ical stationary state with no characteristic time or length scale, without the 
fine tuning of external parameters. We are going to show that the slow driv- 
ing condition is in fact a fine tuning, which make the previous definition 
of SOC ambiguous. A more appropriate definition of SOC, that takes into 
account this observation, is given in p3 |. 



An example of SOC is provided by sandpile models: sand is added grain 
by grain on a rf-dimensional lattice until unstable sand (too large local slope 
of the pile) slides off. In this way the pile reaches a steady-state, in which 
additional sand grains fall off the pile by avalanche events. The steady-state 
is critical since avalanches of any size are observed. This class of models can 
be used to describe a generic avalanche phenomenon, interpreting the sand 
as energy, mechanical stress or heat memory. 

The common characteristic of SOC systems is the presence of a non- 
equilibrium critical steady-state, which we can analyze using the DDRG for- 
malism. We studied the critical height sandpile automaton |R| and the 



forest-fire model [IS, 20]. These two models have been intensively studied 
numerically and can be considered as mile-stones in the field of SOC phe- 
nomena. In what follows we will show that the DDRG allows us to calculate 
analytically the critical exponents and to clarify the SOC nature of both 
models: we are able to study the fixed point and to identify the control and 
the order parameter of the models. 
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4.1 Forest-fire model 



The Forest-Fire model (FFM) has been introduced by Bak et al. |L9| as a an 
example of SOC, and has been then modified by Drossel and Schwabl fl20f . 

The model is defined on a lattice in which each site can be empty (<7j = 0), 
occupied by a green tree (<7j = 1) or by a burning tree (cr^ = 2) (see Fig [l]). 
At each time step the lattice is updated as follows: 

i) a burning tree becomes an empty site; 

ii) a green tree becomes a burning tree if at least one of its neighbors is 
burning; 

iii) a tree can grow at an empty site with probability p\ 

iv) a tree without burning nearest neighbors becomes a burning tree with 
probability /. 

The model was first studied in the case / = for the limit of very slow 
tree growth (p — > 0). In this limit the critical behavior is trivial: the model 
shows spiral-shaped fire fronts separated by a diverging length £ ~ p~ 



where v p ~ 1 |34}j . In the case / > 0, the system was supposed to exhibit 
SOC under the hypothesis of a double separation of time scales: trees grow 
fast compared with the occurrence of lightnings and forest clusters burn down 
much faster than trees grow. This request is expressed by the double limit 
9 = f /p —>■ and p — > 0. The critical state is characterized by a power 
law distribution P(s) = s~ T of the forest clusters of s sites (avalanches in 
the SOC terminology) and the average cluster radius (the correlation length) 
scales as R ~ 9~ u . 

In the past few years, a great amount of work has been done in order 
to describe the critical state of the Forest-Fire model and to calculate the 
critical exponents. Numerical simulations |JR show that in the time scale 
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separation regime the model is close to a critical point with the avalanche 
critical exponent r given by r ~ 1 in d = 1 and r~1.15ind = 2. For the 
exponent u, describing the divergence of the average cluster radius as 9 — > 0, 
it has been found u ~ 1 in d = 1 and v ~ 0.58 in d = 2. The one dimensional 
result has been recovered exactly in [ 37| . Simulations were performed also in 



higher dimensions |3q| . The critical dimension is believed to be d = 6. 

To apply the DDRG to the FFM we follow step by step the strategy 
outlined in Sec. (3). For the sake of simplicity let us first consider in full 
detail the one-dimensional case. To define in practice the DDRG we first 
have to chose a form for the coarse graining operator. We use a cell-to-site 
transformation with scale factor 6 = 2. In this case the operator TZ can be 
written in the following way: 

K{S,a) = l[K(Sj,{<T i }j) (29) 
j 

where each term is acting on a specific cell J and {<Ji}j denotes the configu- 
rations of spins belonging to that cell. Therefore, given a two sites cell, the 
operator TZ renormalizes it in a coarse grained site following only "inside the 
cell" rules. The rules defining TZ £1X6 clS follows. A two sites cell is renormal- 
ized as a tree site if it is spanned from left to right by a connected path of 
green sites. Accordingly, a cell is empty if it is not spanned by a connected 
path of green sites. Finally we consider a cell as burning if it contains at 
least one burning site. In this last case the spanning condition that ensures 
connectivity properties is not necessary because fire spreads automatically 
to nearest neighbor sites. The above renormalization prescription, i.e. the 
operator TZ, is summarized in Fig. |]. We denote with an index a each two 
sites configuration {<Ji}j so that X^o-jj — * J2 a - The dynamical rules of the 
FFM are local, therefore we can define matrix elements reduced to a single 
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sitef] for the dynamical operator T: 

(0 | T | 2) = 1 (30) 

(0 | T | 0) = 1 — p ; (l|T|0)=p (31) 

(1|T|1) = 1-/ ; (2|T|1> = / (32) 

In addition to the above rules, fire spreads between nearest neighbor sites. 
The DDRG second step is the evaluation of the dynamical operator acting 



on the coarse grained variables via the the renormalization equations (|22jj23| ). 
We adopt a finite lattice truncation on the two sites cells subspace defined 
by the operator TZ, and we obtain the single site renormalized dynamical 
operator as 

{Si , rpi , s o } = J2 a J2A^\Tj\a)w a (33) 

where | a) and | a') are the two sites cells statesQ which renormalize respec- 
tively in | Sf) and | Si). We keep the subscript i for the latter states because 
they are referring to a single coarse grained site and not to a system's con- 
figuration. With W a we denote the stationary statistical weight of each a 
configuration. 

Let us now evaluate explicitly the above equations. Because we are in- 
terested in the critical behavior for / < 1 and p < 1 we can write non 
trivial RG equations keeping only terms up to the first order in p and /. In 
addition, in order to define consistently the recursion relations, the operator 
T' must preserve its form at each scale; i.e. no proliferations are allowed. 
This implies that 

(Si = | T' | 5° = 2) = 1 + 0(p\ f 2 ,pf) (34) 

6 1 >, I 1 > and | 2 > are states in which the site i is in the corresponding <n state 
irrespective of the remaining of the system. 
7 For instance a = 1) =| 1, 1) (see fig. 2). 
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where the higher order terms will be neglected in the recursion relations. 
This implies that new parameters are not introduced in the description of 
the system. Thus, to avoid proliferations in the burning event, we have to 
define a dynamical operator T feZ that leaves invariant Eq. fl34|) . It is easy 
to check that this is the case if T b " = T 2 ; i.e. z — 1. In two time steps a 
burning cell evolves always in an empty one (Fig. |3|a) if we neglect second 
order contributions in p and /. In Fig. |||b we show a possible proliferation, 
which however has a weight p 2 , and can therefore be neglected. Avoiding 
proliferations in the burning event defines unambiguously the time scaling 
factor: the relevant time scale in the system is that of the burning process, 



as was already pointed out on the basis of numerical simulations [20] . 

The renormalization of the lightning probability in this framework is 
straightforward. We have only one starting configuration, i.e. a = 1) 
and the recursion relation is 

/' = (51 = 2 | T' I Sj = 1) = E (a' | T 2 | a = 1) = 4/ + 0{p\f\pf) (35) 

q'=4 

In the same way we obtain the expression for the renormalized growth prob- 
ability p as 

from which follows 

p'= 2 n^T^+°(^ 2 ^)- (37) 

A process that contributes to the above equations is shown in Fig. f|a. Finally, 
we have also to treat the normalization of the operator T' . In Fig. |]b we show 
a proliferation given by a process in which with probability / an empty cell 
becomes a burning cell at the coarse grained level. In the previous equations, 
this last process is present only because of the normalization condition: 

(S, = 1 | T> | 5? = 0) = 1 - 2p w ^ w - - /f Jf^ + CV, f.pf) (38) 
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In order to avoid this proliferation we have to impose a supplementary condi- 
tion. If we restrict our analysis to the region in which / <C p, we can neglect 
terms linear in / where terms linear in p are present. This corresponds to 
truncate Eq. (|38|) by keeping only the term linear in p, and thus eliminating 
the proliferation. The inclusion of this proliferation would describe the be- 
havior of the model in the limit f — p which, from numerical simulations, is 
expected to be different. It is worth to remark that the above approximation 
corresponds to renormalize in a separate way the tree growth parameter p 
and the lightning parameter /, assuming that they do not affect each other 
since they act on very different time scales. 

The steady state probability distribution W a is approximate following the 
DDRG general scheme in the lowest order 

W a = {[(p (7l )=n(p (7l )(p (72 ) (39) 
i=i 

where n takes into account the multiplicity due to symmetries of each con- 
figuration. 

We can obtain the densities in the steady state from the following dy- 
namical mean-field equations: 

( Po (t + l)) = (l-p)(p (t)) + (p 2 (t)) (40) 
( Pl (t + 1)) = (p (t))p + (1 - / - (2d - l){p 2 (t))){ Pl (t)) (41) 
(p 2 (t + 1)) = ( Pl (t))(f + (2d - l)(p 2 (t)» (42) 

where d is the spatial dimension (see Ref.|36j for their derivation). The long 

time limit (t — > oo) solution of the above equations^ provides the driving 

condition; i.e. the average steady state densities. 

Collecting all these equations we obtain the DDRG recursion relations 

8 It is worth to remark that also the mean field equations are written for / and p close 
to zero. 
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for the Forest Fire model, that in one dimension read as follows: 



P 
9' 



2p 



(pi) 



(P1>+|(P0>' 
(Pl) + ^(Po). 



(Pl> 

2 



29- 

g(Q,p,p) 
(po) 

(p 2 ) - (1 - ( Pl »a 



(1 - (pi))fl/p = 0; 

a ■ (pi) = 0; 

0. 



(43) 



where we defined a = p/(l +p). We express the recursion relations in terms 
of the parameter 9, in order to compare with numerical simulation. It is 
important to emphasize again that the recursion equations ( |43"D are valid 
only in the double time scale separation f ^ p ^ 1 which defines the range 
of validity for our approximations. This limit is the one for which the FFM 
shows non trivial critical behavior. 

The flow diagram is stable with respect to different coarse graining rules, 
and we find a repulsive fixed point in 9 C = and p c = 0. In order to 
discuss the critical behavior we have to linearize the recursion relations in 
the proximity of this fixed point and to find the relevant eigenvalues of the 
diagonal transformation: 



_ d9' 
Xl ~^9 



8c ,Pc 



dp' 
dp 



(44) 



The fixed point is repulsive, thus defining the critical exponent v in term 
of the largest eigenvalue A of the linearized renormalization equation v = 
log2/log(A) = 1.0. ^From simple scaling relations it is possible to obtain 
also the other critical exponents which are summarized in table I. In this 
respect it is interesting to note that our method yields in the one dimensional 



case the exact results of the rigorous treatment of Ref.|37|]. This is due to 
the relative simplicity of the one dimensional case, where the approximations 
involved - i.e. spanning conditions or proliferations - are irrelevant. 
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In d — 2 the calculation of the RG equations proceeds along the lines 
shown above |TEj but are algebraically more complex. In fact, one has to 
consider the average over different paths, and new dynamical interactions are 
generated at each RG step. This is a signature that we need an approximation 
which truncates the parameter space after each iteration so that it remains 
closed. This is done by considering just the leading order in / and p in the 
renormalization equations, and ignoring any proliferations generated at each 
group iteration. With this scheme we obtain z — 1, which is not an exact 
result also if in good agreement with numerical simulations (z = 1.04 ||38|| ) . 
The fixed point In d = 2 the fixed point remains 6 C = and p c = 0, and 
the largest eigenvalue is given by Ai, which determines the leading scaling 
exponent v = log b/ log Ai = 0.7 (for b = 2). The result is in good agreement 
with numerical simulation (v = 0.6 (081). The numerical value can be further 



improved by using larger cells [18 



It is worth to remark that the DDRG allows to overcome the approxi- 
mations present in the approach of Ref.JT^], where the time scaling was not 
properly considered because of the assumption of an infinite time scale sep- 
aration. In fact, in the limit / = the critical behavior is governed by the 
second eigenvalue A2. This eigenvalue and its relative exponent describes the 
behavior of the correlation length in the deterministic FFM. As opposed to 
Ai, the value of A2 depends on the absolute value of the time scaling factor, 



and therefore could not be obtained by the scheme used in Ref.[18|. The 



numerical value we obtain in d = 1, 2 is u p = log 2/ log A2 = 1.0, which is in 
excellent agreement with the simulation results v v ~ 1 |34| . 



In table I the results obtained in d = 2 for b = 2 are compared with nu- 
merical simulations. The existence of a relevant scaling field and the general 
structure of the flow diagram is stable with respect to different approxima- 
tion schemes, and more refined calculations lead systematically to an im- 
provement in the numerical values of the results. The FFM is critical along 
the line 9 = of the phase space, so that 6 is equivalent to the reduced 
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temperature in thermal phase transitions. In other words 9 is the control pa- 
rameter of the model, and the critical state is reached only by a fine tuning of 
9 to its critical value. The control parameter 9 is the ratio between two very 
different time scales, controlled by / and p, with a critical value fixed to zero. 
In this situation, however, the existence of a time scale separation makes the 
system very close to the critical point without an apparent fine tuning of 
internal parameters. Strictly speaking however, the system is critical just in 
correspondence of the critical value of the control parameter. 

4.2 Sandpile models 

Sandpile models are cellular automata [13|, ^] defined in a d— dimensional 
lattice. A variable E(i), that we denote by energy, is associated with each 
lattice site i. At each time step an input energy SE is added to a randomly 
chosen site. When the energy on a site reaches a threshold value E c the site 
relaxes transferring energy to the neighboring sites: 

E(i)^E(i)-^AE(e) (45) 

e 

E(i + e) -> E(i + e) + AE(e) (46) 

where e represent the unit vectors on the lattice. A typical choice for the 
parameters is , for example, E c = 4 and AE(e) = 5E = 1, but other possibil- 
ities have also been considered. The relaxation of the first site can induce a 
series of relaxations generating an avalanche. Note that the energy is added 
to the system only when the configuration is stable (i.e. all the sites are be- 
low the threshold). The boundary conditions are usually chosen to be open 
so that energy can leave the system. In these conditions the system organizes 
itself in to a stationary state characterized by avalanches of all length scales. 
In particular the distribution for avalanches sizes s decays as a power law 
P(s) ~ s~ r , and the linear size of the avalanche scales with time r ~ t z . 
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This model has been extensively studied in the past by means of numerical 
simulations [HDL |4"I|, [4"2"[ and several exact results have been derived for Abelian 



sandpiles 

In the steady state each configuration of the system can be described by 
the energy probability distribution W({Ei}). In the zero order approxima- 
tion we consider the incoherent part of the distribution: 

^ (i) ({^})=n^(^) ( 47 ) 

i 

where w(E) is the single site energy distribution. The average energy of a site 
evolves according to the following equation, written in continuum notation: 

^ = SE m - SE out (48) 

where 5Ei n is the average energy entering into the site either because of a 
relaxation in a neighboring site or because of the external perturbations, and 
5E out is the average energy dissipated by the site. At lower scale 
SE out = (p) ^ e AE(e), where we have defined 

(p) = w{E)dE (49) 



IE C -5E 

as the probability that a site relaxes. 

This equation suggest a simple way to describe the steady state of the 
model. At any scale, we can divide the sites in critical (<7j = 1) and stable 
((7j = 0). Stable sites do not relax when energy is added to them. On the 
other hand critical sites relax when they receive an energy grain SE in . In this 
formalism p represent the density of critical sites. For convenience we will 
also define unstable sites {&{ = 2), as those that are relaxing, even though 
they are not present in the static configurations of the system (see Fig. [I]). 
These definitions can be extended to a generic scale b. For instance, a cell 
at scale b is considered critical if the addition of energy 5E in {b) induces a 
relaxation of the size of the cell (i.e. the avalanche spans the cell). 
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The DDRG allows us to develop a general renormalization scheme for 
sandpiles which put in a broader and more systematic context the approach 
of Ref . [jr7| . To construct the DDRG recursion relations we have to describe 
the relaxation at a coarse grained scale. We first note that the only non 
trivial matrix element of the time evolution operator is the one describing 
the evolution of a cell from unstable to stable. This process occurs when 
a critical cell becomes unstable due to the addition of energy. The inverse 
process, a stable cell becoming critical, is simply due to the accumulation 
of energy and is not characterized by critical exponents which, on the other 
hand, describe the avalanche propagation. 

In a relaxation at the minimal scale energy is distributed equally in the 
four directions. This is no longer the case at a coarse grained level where 
different possibilities arise: the energy in principle can be distributed to one, 
two, three or four neighbors. It is also worth to remark that in certain 
case unstable sites at the coarse grained scale do not dissipate energy to 
nearest neighbors, representing just intra-site energy rearrangements. These 
processes defines the probability that relaxation events take place on the 
renormalized scale without energy transfer. All these events occur with prob- 
abilities 

P = (Po,Pi,P2,P3,Pa) 

In terms of the matrix element (0|T|2) the vector P represents the probabil- 
ities 

Pn = (0|T|2>„ (50) 

where (0|T|2) n is the probability that a relaxing site becomes stable and 
transfers energy to n neighbors. In this way we have obtained the set of 
parameters that describes the dynamics. Of course the choice of the parame- 
ters space is not uniquely determined, one encounters proliferation problems 
typical of real space RG methods. For instance, higher orders proliferations 
are due to multiple relaxations of the same site and sites becoming critical 



Vespignani, Zapperi and Loreto 



29 



during the dynamical process (i.e.: (1|T|2)). In the following the practical 
implementation of the method considers just the minimal proliferation we 
reported above. 

The renormalized matrix element is then obtained by considering all the 
processes that span the cell and transfer energy outside. This rule implicitly 
defines the effective dynamical operator T fez (see App.A), the renormalized 
parameters being: 



We proceed in defining explicitly a renormalization procedure for the dynam- 
ics by considering a finite truncation on four-sites cells. This corresponds to 
a cell-to-site transformation on a square lattice, in which each cells at the 
coarser scale is formed by four sub-cells at the finer scale: the length scaling 
factor is b = 2. The operator 1Z which define the coarse grained variables 
acts on each specific cell through "inside cell" rules. A cell is renormalized 
as a relaxing one if it contains a relaxing sub-cell which transfers energy to a 
critical sub-cell. In this way we ensure that the occurring relaxation process 
is extending over the size of the renormalized length scale independently of 
the successive avalanche evolution. A critical cell is therefore defined by a 
cell which can be spanned by a path of relaxation events. The scheme con- 
siders only connected paths that span the cell from left to right or top to 
bottom. This spanning rule implies that only paths extending over the size 
of the resulting length scale contribute to the renormalized dynamics, and it 
ensures the connectivity properties of the avalanche in the renormalization 
procedure. 

Every cell at the coarser scale can be characterized by an index a that 
indicates the configuration of sub-cells, and we have that ~> J2 a - The 

approximated stationary distribution (Eq. fl2lD ) for each of these configura- 



Pn = 



E W Q } E W} nS° = 2,a°)^(g = 0,a)(a | T» z | a°) n W(a ) 



(51) 
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tions is given by: 

W a ((p)) = n a Y[(p ai ) (52) 

i=l 

where n a is a factor due to the multiplicity of each configuration. 

By using this scheme and replacing sums over configurations with sums 
over the index a, the recursion relations can then be rewritten in the simpler 
form 

Pn = ^E^((P))E(«' I T " Z GV) I «)n (53) 

JV a a' 

where | a), \ a') denotes the four sites configurations which renormalize in 
| Sf = 2) and | Si = 0), respectively. In the above expression the denominator 
of eq. ( |5TD is adsorbed in the normalization factor A/". 

The time scaling factor and the explicit definition of the effective dynam- 
ical operator T bz can be found in the appendix A. It is worth to remark that 
the present case is very different from the FFM and the definition of the 
effective dynamical operator is non trivial. The driving condition (eq.(p3|)) 
is obtained from eq.(|58|) by imposing stationarity. This implies that the sta- 
tionary state is characterized by the balance between the energy that goes 
in and the energy that goes out of the system. We assume that energy is 
transferred in "quanta" SE = SE in in each direction and we obtain on average 

6E=(p)6Ej2n Pn (54) 

n 

which implies 

(P) = (55) 

This relation gives the average density of critical sites in the steady-state, 
allowing us to evaluate the approximate stationary distribution at each scale. 
Therefore, eq.s(|53D and (|55|) are the complete set of the DDRG recursion 
relations. The practical calculation of all the paths involved in the evaluation 
of the above equation is very laborious and can be found elsewhere |44|, [4^ . 
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moreover, we have shown a particular truncation scheme and more generally 
the explicit evaluation of the recursion relations depends upon the chosen 
spanning condition and number of proliferations considered. 



In Ref. ||17|| , it has been developed the simplest closed renormalization 
scheme which neglects in addition to the higher order proliferations men- 
tioned above also the probability po. The flow diagram shows an attrac- 
tive fixed point: the parameters evolves spontaneously towards their critical 
value. From the fixed point the critical exponents can be computed. For a 
square cell of size 6 = 2 the results obtained are r = 1.25 and z = 1.17, in 



good agreement with computer simulations [40, 41]. The same method has 



also been applied to dissipative sandpile models |17]] and to directed sand- 



pile models fHfl . The effect of the po processes is being included and the 



results will be reported in a forthcoming paper. Recently, the expressions 
for the recursion relations have been linked to a branching mechanism that 
allows their calculation through a generating function. Using this method, 
it is possible to include more proliferations to the set of relaxation processes 
considered. The results obtained with this improved scheme allows an excel- 
lent qualitative and quantitative description of critical height sandpile model 

It is worth to remark that in this scheme, as is usually done in com- 
puter simulations, we are implicitly assuming a slow driving condition for 
the model. In fact, in the evaluation of the RG equations the external drive 
does not interfere with the relaxation process: dynamical processes during 
which sites become critical are not considered. In order to overcome this 
approximation the energy flow on its turn should be renormalized and the 
addition of energy (possibility of new relaxation events) during relaxation 
processes be allowed. Enlarging the phase space, a relevant parameter ap- 
pears, e.g. the driving rate. The driving rate is the incoming current of 
energy per unit time with respect to the total average energy in the system. 
Also in this case, as in the FFM, the system would be strictly critical just 
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in correspondence of infinitesimal driving rate. ^From this point of view the 
only difference between Sandpile models and Forest-Fire models is that the 
latter can not be studied in a subspace with no relevant parameters without 
destroying the model itself. The meaning of the self-organization is then re- 
lated to the widespread existence of systems with very different time scales 
and not to the absence of relevant control parameters as often reported in 
literature. 

5 DDRG and equilibrium critical phenomena 

The DDRG represents a general method to approach non-equilibrium critical 
systems with a stationary state and it allows also to study equilibrium models 
at the critical point. In this last case the stationary state is characterized by 
probability densities written in terms of the Hamiltonian of the system by 
means of the Gibbs distributions. This kind of systems have been extensively 
studied and well established theoretical tools are available to approach them. 
The prototype of such systems is the Ising model. On each site % of some 
finite dimensional lattice we place a random variable <7j taking the values ±1. 
The Hamiltonian is 

H(v) = -Jj2^ a j ( 56 ) 

where the sum runs over all the nearest neighbors pairs. The correspondent 
Gibbs measure is 

W(<t) = Z~ x exp {-(3H(a)} (57) 

where Z is a normalization factor (the partition function). 

In order to apply the DDRG to equilibrium critical phenomena one has 
to treat these systems as dynamical systems with a well-defined stationary 
state. Our discussion will refer to the Ising model with heat bath dynamics 
but the arguments can be generalized to different cases. In a system described 
by a state {<7j}, the conditional probability to flip a spin in a site % is given 
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by 

where nn indicates the nearest-neighbors of the site i and (3 is the inverse of 
the temperature. The configuration {o^} is obtained from {of} by flipping 
the spin %. The heat-bath algorithm updates of by choosing a new spin value 
<jj, independently of the old value of of, from the conditional probability 
given by Eq. (|58|). All the others spins remain unchanged. If {o^} and {of} 
are two arbitrary configurations of the system and (a\T(/3)\a°) , (o°|T(/3)|o) 
are the 1-step transition probabilities between the two configuration, the 
following relation, known as detailed balance, is satisfied by the heat-bath 
dynamics 

W(a)(a\T(P)W°) = W(o°)(o°\T(l3)\v) ( 59 ) 

where W(a) and W(a°) are the stationary (equilibrium) distributions for the 
states {oj} and {of}. 

^From the detailed balance we can deduce a stationarity condition for 
the dynamics which provides the steady-state distribution to be used in the 
DDRG. For the actual implementation of the RG procedure one can use 
standard techniques for the definition of the operator 7Z, e.g. majority rules, 
spanning conditions etc., and for the dynamical renormalization of the pa- 
rameters [i. The work in this direction is still in progress and it will provide 
a test for the flexibility of our approach. In order to test the accuracy of the 
method we could compare the results for the probability densities obtained 
in our framework with those obtained renormalizing directly the Gibbs dis- 
tributions with usual RG schemes. 



6 Conclusions. 



In this paper we have presented a new renormalization scheme especially 
suited for systems with non-equilibrium critical steady-state. The essential 
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idea of the method is the use of an approximate stationary probability dis- 
tribution for the configurations of the system. This distribution is evaluated 
through a driving condition which identifies the single time averages in the 
steady-state. This approximate distribution is used as a weight in the renor- 
malization of the master equation which takes into account correlations due 
to the dynamical evolution. These correlation are thus considered in the 
approximate probability distribution of the renormalized system, which is 
calculated by the driving condition with renormalized dynamical parame- 
ters. The dynamical renormalization of the master equation is based on real 
space RG schemes which, in spite of their low systematicity, are very sim- 
ple and intuitive. In addition, these schemes leave room for "higher order" 
techniques to improve the accuracy of the results. 

The DDRG can be used to study systems with stationary critical state, 
but is particularly useful in non-equilibrium systems for which it is impossible 
to obtain the stationary probability distribution in the configuration space. 
In fact, the DDRG appears to be extremely powerful in the study of SOC 
systems. These systems show a non-equilibrium steady-state very close to 
the critical point for a wide range of internal parameters. 

The application to forest-fire and sandpile models shows that the gen- 
eral results obtained are not affected by the approximations involved in our 
scheme, even though a more refined treatment of the calculation scheme, i.e. 
by introducing more proliferations, leads to an improvement of the values 
obtained for the critical exponents. For SOC models the DDRG can be con- 
sidered as a general theoretical framework that provides the analytical tools 
for a qualitative and quantitative study of the critical stationary state. 

We considered as well the application to usual equilibrium phenomena, 
in the perspective of a better understanding of the method and its applica- 
tion to other systems. In particular the DDRG could be very effective to 
study the critical properties of driven diffusive systems || [TO], which to our 
knowledge have never been approached by real space RG methods. Work in 
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this direction is currently in progress. 
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Appendix A 

The explicit evaluation of the time scaling factor and the corresponding 
recursion relation is a complex task, because there is no standard recipe 
for this in real space RG schemes. Here we show a general strategy, which 
depends on the specific model for its actual implementation. 

The effective operator T h " contains all the dynamical processes that con- 
tribute to the definition of a meaningful renormalized dynamics. We define 
the following transformation 

(*\T b \ri\a ) = Y,V N (a\T N (ri\a ) (A. 1) 

N 

where V N is the renormalization operator for the dynamical evolution of the 
system: it is a projection operator that samples only the paths of A" time 
steps which have to be considered in the definition of the effective operator 
T bZ . To clarify this point, let us consider for a moment a spin flip dynamics 
in a Ising-like system. The matrix element T(/x) is non zero only for those 
configurations that are related by a single spin flip. In order to preserve the 
same form for T at a coarse grained scale, we have to impose that the renor- 
malized time evolution operator connects only configurations that differ by a 
single coarse grained spin flip. Flipping a coarse grained spin corresponds to 
the subsequent flipping of different spins in the original system. The number 
of flipping necessary to flip a macroscopic spin is not uniquely defined, but 
depend on the configuration, both at the coarse grained and fine grained 
level, and on the particular dynamical path chosen. The effective operator 
T bZ is then a convolution of different A^-steps operators. 

The operator T>m is chosen on the basis of physical considerations: span- 
ning conditions etc. In addition, TV should satisfy some general properties 
in order to preserve the symmetry or the internal space of the dynamical 
variables. For instance, we have to ensure the normalization of the effective 
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dynamical operator by the property 

J2J2V N (a\T N (»)\a°) = l (A. 2) 

M n 

Moreover, T>n must be consistent with the definition of the renormalization 
operator 1Z: it should describe dynamical processes among renormalized vari- 
ables of the same type of those given by the operator 1Z. Finally, T>n has to 
preserve the form of the dynamical operator T at each scale. This condition 
imposes that the time scaling is consistent with the length scaling used in 
1Z. In this way it is possible to map the renormalized system in the old one 
with renormalized variables. 

As previously mentioned the operator T>n can assume a very simple form: 

£V = $n,n' (A. 3) 

where N' = b z . In general, however, more complicated expressions are en- 



countered (see Ref. [[17], |18|), since T>n depends on the specific dynamics. We 
have defined the effective evolution operator so that T'(/i) = T(//) for the 
renormalized system. On the other hand, the operator T b * is in general the 
convolution of the discrete time step operator T N projected by the renormal- 
ization operator T>^. Thus, we can write it as a sum of terms that represent 
the statistical weight for the evolution paths of N time steps: 

W\T bZ (v)\a ) = j:^Ari ( A - 4 ) 

N 

These terms, which are obtained by the specific definition of the operator 
V N , can be used to calculate the time scaling factor b z = g(u) as an average 
over the renormalized dynamical processes. In fact, we can interpret the right 
sum in Eq. ([A. 4 ) as an integral over dynamical evolution paths of different 



time duration; the time scaling factor being an average over these paths 
whose statistical weights are the terms 7^0 (a*)- For each case we have to 
find the process which defines the relevant time scale of the phenomenon. 
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The time scaling factor will be the average over contributing paths to the 
renormalization of this process. 

To illustrate how the above procedure works in practice, we shall discuss 
explicitly the FFM and sandpile examples. For FFM we are in the simple 
case of Eq. flA. 3|) . In fact the relevant process, i.e. the burning process, 
determines Vjsr = 8^,2- therefore only terms J^ a oip, /) = (er | T 2 (p,f) \ a ) 
are allowed, so we simply have b z = g(p, /) = 2. 

For the sandpile automata the calculation is rather laborious. In the slow 
driving regime, the only relevant time scaling length is given by the single 
relaxation event, for which we have at each scale 

(Si = | T' | S? = 2) = 1 (A. 5) 

The operator T>^ selects those spanning paths that lead to a relaxation pro- 
cess at the coarser scale. By using the cell-to site transformation defined in 
Sec. 4. 2, it is easy to show 

Y,K(S° = 2,a°)K(S i = 0,<7)£ 7 * o(?v) = EE^W I T*M I «> 

{o-} N a' N 

(A. 6) 

where | a), \ a') denotes the four sites configurations which renormalize in 
| Sf = 2) and | Si — 0), respectively. The sum over n denotes we are con- 
sidering relaxation processes without distinguishing the number of affected 
nearest neighbors. The operator T>^ is therefore defined explicitly as an op- 
erator acting on the paths internal to four sites cells. It selects for each 
just relaxation paths which consist of A^ connected non-contemporary relax- 
ation events that leave the cell without unstable sites. In a mathematical 
forms it reads as 

JV-l 4 

V N = I] (1 " <VJ II E S(m - E W ( A - 7) 

ie{a'} J=0 m=l iefaj} 

where aj's are the intermediate cell configurations during the dynamical 
evolution and J2ie{ aj } denotes the sum over all the sites in the cells. In the 
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above expression, each delta function acts on a different intermediate cell 
eliminating those paths which do not have activity at each dynamical step. 
Furthermore, the operator ensures that in the cell a' (Nth step) no activity 
is present; i.e the process has stopped. Finally we have to write the equation 
that gives the time scaling factor from the total average over contributing 
processes to the renormalized matrix element (0 | T' | 2) : 



where we used the DDRG scheme to explicitly get the stationary weights and 
M is an opportune normalization factor. 

The above relations will provide the consistent rescaling of time by im- 
posing that b z = g(Pn) from which it is possible to calculate the dynamical 
critical exponent. This also shows that in general the factor g(fi) is a function 
of the dynamical parameters. 



#0) 



Z { aO } 1Z(S? = 2,a°)W(a°) 



(A. 8) 



and by inserting Eq.( |A. 6| ) in the above relation we finally obtain 




(A. 9) 
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FIGURE CAPTIONS 

• Fig.l The sites classification in the Forest Fire and in the sandpile 
model. 

• Fig. 2 The rules defining the renormalization operator 1Z for the one 
dimensional Forest Fire model. 

• Fig. 3 (a) Evolution of a burning cell into a empty cell in two time 
steps, (b) A possible proliferation in which a burning cell evolves in 
two time steps into a green cell. This process is of order p 2 and can 
therefore be neglected. 

• Fig. 4 (a) A process contributing to renormalization of p: an empty cell 
becomes green due to the growth of one site, (b) A possible proliferation 
of the growth dynamics in which an empty cell becomes burning. In 
the limit / <C p this process can be neglected. 
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TABLE CAPTIONS 

• Table I In this table we summarize our results for the Forest Fire Model 
critical exponents obtained with different approximation schemes. For 
comparison we report also the exact or numerical results. *: Exact 



results from |p7[| . . Numerical results from |35|, ^8 . 
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d = 1 
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RG 


1.0 


1.0 


1.0 


Exact results* 


1.0 


1.0 


1.0 


d = 2 


V 


z 
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RG 2x2 


0.73 


1.0 


1.19 


Numerical results" 4 " 


0.58 


1.04 


1.15 



TABLE I 
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